An introduction to quantum annealing 

Diego de Falco and Dario Tamascelli 

Dipartimento di Scienze dellTnformazione, Universita degli Studi di Milano 

Via Comelico, 39/41, 20135 Milano- Italy 

e-mail:defalco@dsi. unimi.it, tamascelli@dsi.unimi.it 



in 



cr 



> 



c3 



Abstract 



Quantum Annealing, or Quantum Stochastic Optimization, is a 
classical randomized algorithm which provides good heuristics for the 
solution of hard optimization problems. The algorithm, suggested by 
the behaviour of quantum systems, is an example of proficuous cross 
contamination between classical and quantum computer science. In 
this survey paper we illustrate how hard combinatorial problems are 
tackled by quantum computation and present some examples of the 
heuristics provided by Quantum Annealing. We also present prelimi- 
nary results about the application of quantum dissipation (as an alter- 
native to Imaginary Time Evolution) to the task of driving a quantum 
^. . system toward its state of lowest energy. 

q 

q ■ 1 Introduction 

Quantum computation stems from a simple observation [23j H2! : an y com- 
putation is ultimately performed on a physical device; to any input-output 
relation there must correspond a change in the state of the device. If the 
5_i ■ device is microscopic, then its evolution is ruled by the laws of quantum 



mechanics. 

The question of whether the change of evolution rules can lead to a break- 
through in computational complexity theory is still unanswered |11 [ I34 [ |40"]. 
What is known is that a quantum computer would outperform a classi- 
cal one in specific tasks such as integer factorization [47J and searching an 
unordered database |28[ 129], In fact, Grover's algorithm can search a key- 
word quadratically faster than any classical algorithm, whereas in the case of 
Shor's factorization the speedup is exponential. However, since factorization 



is in NF^J but not believed to be NP-complete, the success of factorization 
does not extend to the whole NP class. 

Nonetheless, quantum computation is nowadays an active research field and 
represents a fascinating example of two way communication between com- 
puter science and quantum physics. On one side, computer science used 
quantum mechanics to define a new computational model, on the other the 
language of information and complexity theory has allowed a sharper un- 
derstanding of some aspects of quantum mechanics. 

Moreover, it is clear that quantum mechanics will sooner or later force its 
way into our classical computing devices: Intel® has recently (Feb 201Co) 
announced the first 25nm NAND logical gate: 500 Bohr radii. 

In this paper we present an introduction to quantum adiabatic computation 
and quantum annealing, or quantum stochastic optimization, representing, 
respectively, a quantum and a quantum-inspired optimization algorithm. In 
Section [2] we introduce the class of problems that are hard for a quantum 
computer. We consider the problem 3-SAT, its mapping into a quantum 
problem and the quantum adiabatic technique. In Section [3] we describe 
Quantum Annealing, a classical algorithm which captures the features of 
the quantum ground state process to provide good heuristics for combinato- 
rial problems. Section [4] is devoted to dissipative quantum annealing and to 
the presentation of preliminary heuristic results on the dynamics of a dissi- 
pative system. Concluding remarks and possible lines of future research are 
presented in the last section. 

2 Transition to the quantum world 

The need of considering the features of quantum computing devices led to the 
definition of the quantum counterpart of some classical complexity classes. 
Quantum algorithms are intrinsically probabilistic, since the result of the 
measurement of an observable of a quantum system is a random variable. 
Not surprisingly, the quantum correspondent of P is BQP (Bounded-error 
Quantum Polynomial-time), the class of decision problems solvable in poly- 
nomial time by a quantum Turing machine, with error probability at most 
1/3. 

QMA (Quantum Merlin Arthur) is the quantum counterpart of the class 



lr To be more precise, factorization is in FNP, the function problem extension of NP. 
The decision problem version of factorization (given an integer N and an integer M with 
1 < M < N, does N have a factor d with 1 < d < Ml) is in NP. 

2 http://www. intel.com/pressroom/archive/releases/2010/20100201comp. htm 



NP|50|. It is defined as the class of decision problems such that a "yes" 
answer can be verified by a 1-message quantum interactive proof. That is: 
a quantum state (the "proof") is given to a BQP (i.e. quantum polynomial- 
time) verifier. We require that if the answer to the decision problem is "yes" 
then there exists a state such that the verifier accepts with probability at 
least 2/3; if the answer is "no" then for all states the verifier rejects with 
probability at least 2/3. 

All the known complete problems for the class QMA are promise problems, 
i.e. decision problems where the input is promised to belong to a subset 
of all possible inputs. An example of QMA-complete problem is the k-local 
Hamiltonian problem: 

INSTANCE: a collection {Hi,H2, ■ ■ ■ ,H n } of Hamiltonians each of which 
acts on at most k qubits; real numbers a, b such that b — a = 0(l/poly(n)). 

QUESTION: is the smallest eigenvalue of Y21=i Hj ^ ess than a or greater 
than 6, promised that this is the case? 

It has been shown that this decision problem is complete for the class QMA 

even in the case of 2-local Hamiltonians [34 1. 

The well known problems A;- SAT, which are complete for NP for k > 3, can 

be encoded in a fe-local Hamiltonian problem. 

For instance, let us consider a given boolean formula: 

C\ A C 2 A . . . A C M 

over the variables x±, X2, • • • , xn] each C{ is the disjunction (b^ V b-i 2 V . . . V bi h ) 
and bi is either Xi or its negation -ixj. Our problem is to decide whether 
there is an assignment to the boolean variables x\,X2, ■ ■ ■ ,xn that satisfies 
all the clauses simultaneously. 
Given a clause Q = (6j a V 6j 2 V . . . V bi k ) we define the local function: 



fc i (x 1 ,x 2 ,...,x N )= JJ(l-a^)- ]~| 



■ ' j 



A!_, A^ containing all the indices of the variables that compare negated or 
not negated in the clause C{ respectively. Given an assignment (x\, x 2 , ■ ■ ■ , xn), 
the function fc i will return if the clause Cj is satisfied, 1 otherwise. 
The cost V(xi,X2, ■ ■ ■ ,xjy) of an assignment can be defined therefore in 
terms of the functions fc t as: 

M 

V(x 1 ,x 2 ,...,x N ) = ^2fOi(xi,X2,...,x N ), 

i=l 



that is the total number of violated clauses. An instance of fc-SAT is there- 
fore satisfiable if there exists a configuration of zero cost. 

The translation of a fc-SAT problem to the corresponding fc-local Hamilto- 
nian, problem is straightforward. 

To each variable rr, we assign a qubit, i.e. a 2-level quantum system. For the 
sake of definiteness we will consider spins (<Ji(i), 02(2), 03(1)), i = 1, 2, . . . , N 
and their component 03 (z) along the z axis of a given reference frame as com- 
putational direction. Furthermore, we decide that if an assignment assigns 
the value 1 (true) to the variable Xi then the corresponding qubit is in the 
state crs(i) = +1 (spin up) and in the state 03 (i) = —1 otherwise. 
To each clause we associate the Hamiltonian term: 

e,. n ^f^- n 1 - ± j M - m 

which acts non trivially only on k qubits. The total Hamiltonian: 

M 
H = Y,Hi (2) 

is fc-local and plays the role of a cost function: if an assignment satisfies all 
the clauses simultaneously then the corresponding energy is zero. Otherwise 
the energy would be equal to the number of violated clauses, thus > 0. 

The fc-local Hamiltonian problem associated to /c-SAT is a particular case 
of the &-QSAT problem, which is formulated as follows: 

INSTANCE: Hamiltonian H = Ysf=\ p j acting on N qubits, where each Pj 
is projector on a 2 k dimensional subspace of the whole (2 ) dimensional 
Hilbert space. 

QUESTION: Is the ground state energy Eq of H zero, promised that either 
E = or E > l/poly(N)l 

If each Pj projects on an element of the computational basis, we obtain 
the fe-local Hamiltonian associated to fe-SAT by the construction described 
above. 

Interestingly enough it is proved that 4-QSAT is QMA complete whereas 
2-QSAT is in P|14|. For 3-QSAT the answer is not known. 

Quantum Adiabatic Computation |21| was welcomed by the quantum com- 
puting community due to the preliminary good results produced when ap- 
plied to small random instances of NP-complete problems |21|. |22| I5T] and 



its equivalence to the quantum circuital model pQ. 

The computational paradigm is based on the well known adiabatic theorem 
[1'6\ [5] . Given a time T > and two Hamiltonians Hj and H? we consider 
the time dependent Hamiltonian: 

H(t) = tH T + (T - t)Hi, 0<t<T, (3) 

or, equivalently, 

H(s) = h(±Y 0< s <1. 

Let us indicate with | s; efc(s) ) the instantaneous eigenvector of -ff(s) corre- 
sponding to the instantaneous eigenvalue ej.(s) with eo(s) < ei(s) < . . . < 
e n (s), < s < 1 and by | ^»(s) } the solution of the Cauchy problem: 

\|V(0)) = |0;e (0)}. 
The adiabatic theorem states that, if the time T satisfies 

T %>-£-, (4) 

where g m in is the minimum gap 

g min = min (ei(s) - e (s)) 

0<s<l 



and 



max 

0<s<l 



s;ei{s) \ — \ s;e {s) 



then |( 1; eo(l) | ipiX) )\ can De m ade arbitrarily close to 1. In other words if 
we start in the state 1 0; eo(0) } we will end up in the ground state | T; eo(T) } 
of the target Hamiltonian Ht- In practical cases £ is not too large; thus 
the size of T is governed by g~ in : the smaller g m in the slower must be the 
change rate of the Hamiltonian if we want to avoid transitions (the so called 
Landau-Zener transitions [53J) from the ground state to excited states. 

The adiabatic method can be used to find the unknown ground state of an 
Hamiltonian Ht- We can start from the known (or easy to prepare) ground 
state of an auxiliary Hamiltonian Hj and consider the time dependent con- 
vex combination ([3]). 



For example let us consider the 3-SAT problem |22j: we need to understand 
if the ground state of the Hamiltonian ([2]) has energy or not. We set 
the target Hamiltonian H? to H. As initial Hamiltonian we consider the 
Hamiltonian 

i 

having a ground state easy to prepare (all the spins aligned along positive 
x direction). 

The time T required to fulfill the hypotheses of the adiabatic theorem is 
then the cost in time of the algorithm and is determined by the gap g m in- 

The results obtained in the seminal paper by Farhi et al. [22J on small 
random instances of 3-SAT suggested that the gap g m in scaled polynomially 
with the problem size. However, subsequent results [Ml SU [Ml E] proved 
that g m i n can be exponentially small. It has been even shown that adiabatic 
quantum computation can perform worse than other heuristic classical and 
quantum algorithms on some instance of 3-SAT |30j. 

Though these results do not prove that 3-SAT is QMA-complete, they hint 
nevertheless that 3-SAT is hard for quantum computers as well. Recent 
studies focused on determining why 3-QSAT has instances which are hard 
for adiabatic quantum computation and what general features they possess. 
The idea is to use typical, but not necessarily worst, random cases with 
a specific density of clauses a = -^ (M number of clauses, N number of 
variables) and to investigate how g m in varies (on the average) as a function 
of a. We refer the interested reader to [40 [ [39j [6] and references therein. 

3 Quantum annealing 

So far we discussed if and how quantum computational devices could be 
used to tackle problems which are hard to solve by classical means. In this 
section we present a class of classical heuristics (i.e. algorithms meant to 
be run on classical computational devices) suggested by the behaviour of 
quantum systems. 



Many well-known heuristic optimization techniques [43] are based on natu- 
ral metaphors: genetic algorithm, particle swarm optimization, ant-colony 
algorithms, simulated annealing and taboo search. In simulated annealing 
[35| . for example, the space of admissible solutions to a given optimization 
problem is visited by a temperature dependent random walk. The cost func- 
tion defines the potential energy profile of the solution space and thermal 



fluctuations avoid that the exploration gets stuck in a local minimum. An 
opportunely scheduled temperature lowering (annealing), then, stabilizes 
the walk around a, hopefully global, minimum of the potential profile. 

The idea of using quantum, instead of thermal, jumps to explore the solution 
space of a given optimization problem was proposed in Refs.|9, 8p. It was 
suggested by the behaviour of the stochastic process [2j [20] q v associated 
with the ground state (state of minimal energy) of a Hamiltonian of the 
form: 



H v 



,2 Q2 



+ V(x), 



(5) 



2 dx 2 
where the potential function V encodes the cost function to be minimized. 




Figure 1: Thermal jumps, which in SA allow the exploration of the solution 
space, are substituted, in QA, by quantum jumps (tunneling). 

For any fixed v, once given the ground state ip u of H u , the stochastic process 
q u can be built by the ground state transformation [2]: let ip„ € L 2 (R,dx) 
be the ground state of the Hamiltonian ([5]). Under quite general hypotheses 
on the potential V, the ground state ip v can be taken strictly positive. The 
transformation: 



U :^^XJi) 



i>v 



For the connection with the QAC of the previous section we refer the reader to 



or ground state transformation, is well defined and unitary from L 2 (R,dx) 
to L 2 (R,ip 2 dx). Under this transformation, H u takes the form: 

H u = UHvU- 1 = -vL v + E„ 

where 

r 1 d2 u 

U = 2 U dx^ + K 

has the form of the generator of a diffusion process q u on the real line, with 
drift 

The behaviour of the sample paths of the stationary ground state process 
q u is characterized by long sojourns around the stable configurations, i.e. 
minima of V(x), interrupted by rare large fluctuations which carry q v from 
one minimum to another: q v is thus allowed to "tunnel" away from local 
minima to the global minimum of V(x) (see Figj2]). The diffusive behaviour 
of q v is determined by the Laplacian term in (j^J), i.e. the kinetic energy, 
which is controlled by the parameter v. The deep analysis of the semi- 
classical limit performed in [32] shows, indeed, that as v \ "the process 
will behave much like a Markov chain whose state space is discrete and given 
by the stable configurations". 

Quantum annealing, or Quantum Stochastic Optimization (QSO), in the 
original proposals of Refs.[9l |8], did not intend to reproduce the dynamics 
of a quantum mechanical system (a task computationally untractable |23|). 
but rather to simulate the ground state process of the Hamiltonian H u as 
!/\0. The algorithm is completely classical and intends to capture the 
features of the quantum process which can allow an efficient exploration of 
the solution space described above. 

The desired ground state estimation is obtained by means of imaginary time 
evolution, that is by letting an arbitrarily chosen normalized initial state 
tp(x, 0) of the system evolve not under the action of e~ " but under e~ ". 
This replacement of the Schrodinger equation by a heat equation has the 
following property: 

li m _ e -*f i |^(o)) = \ 1 j )v ) ) 
t— >oo at 

where at = { 4>o I V'(O) ) exp(— tE'o). The proof of this fact is straightforward: 
given the eigenstates <j>o(= ip u ),<pi, . . . ,4>n of H u and the corresponding 




Figure 2: A sample path of the ground state process q v for small v. The po- 
tential profile is reproduced on the vertical axes. In the inset we reproduce 
the potential V(x) (dashed line) and the ground state for the this example 
(solid line). For the example chosen the probability distribution is more 
concentrated around a local rather than the global minimum. What mat- 
ters, however, is that the process jumps between stable configurations, the 
absolute minimum included. 



eigenvalues Eq < E± < . . . < En, we have: 

lim — e-'^-|V(0)) " ' 

t— >oo at 



lim - , 

t^oo ( O I ^(0) ) 



N _ t 2lL 

i-A-i^X^l^o)) (6) 

n=o e » 

N 



£ 



t E„-E 



lim / Yl e"*^^ 1 1 <t>n ) ( <f>n I ^(0) } 

t^oo ( Q I -0(0) ) ^ 



n=0 



N 



hm | 0o ) + , , ,* E e"'^ I ^ )( K | V-(O) )• 

t->oo ( 0o | ^(0) ) ^ 



The excited states are "projected out" and, asymptotically, only the ground 
state survives. 

We are nevertheless left with the problem of determining e~~ ". A direct 
diagonalization of H v is impossible, given the dimensionality of the operator 
which scales exponentially with the dimension of the system. However, the 
Feynman-Kac formula |33^ [23] : 

{e~ tHv if) (x) = E ( exp (- I F(e(r))dr J if)(e(t)) e(0) = x\ , (7) 

suggests a way to estimate the evolved state e~ "i/> for any initial state 
ip: the component x of the state vector corresponds to the expected value 
of (the exponential) of the integral of the potential (cost) function V along 
the stochastic trajectories e(r). Each trajectory e(r) starts at x and, in 
the time interval (0,£), makes N(t) transitions toward nearest-neighboburs 
with uniform probability; the number of random transition N(t) is a Pois- 
son process of intensity v. A detailed description of the sampling and of 
the ground state estimation process is beyond the scope of this paper and 
we refer the interested reader to Ref.[8J. Here it suffices to say that once 
given an estimate of ip u {y) for every neighbour y of the current solution x, 
the exploration proceeds with high probability toward the nearest-neighbour 
solution having the highest estimated value of ip u . In Quantum Stochastic 
Optimization, therefore, each move is local (i.e. to a nearest-neighbour of 
the current solution) but the decision rule on which neighbour to accept is 
based on the prospection of an ensemble of long chains e(r). 
When implemented in a working computer program, the procedure described 
above requires a variety of approximations. First of all, equation ([7]) repro- 
duces the ground-state only in the limit t — > oo, corresponding to sample 
paths of infinite length. The expected value of the right hand side of ([7]), 
moreover, will be estimated by means of a finite size sample. Finally, the 
number of neighbours of a configuration may be too large to allow the es- 
timation of i/j v on the whole neighbourhood of a given solution. The accu- 
racy of the approximation depends, for example, on the actual length vt of 
the sample paths, the number of paths n per neighbour and the dimension 
| Neigh\ of the subset of the set of neighbours of a given solution to consider 
at each move. Some "engineering" is also in order [9] (see the pseudo-code of 
Quantum Stochastic Optimization reported here): for example, we can call 
a local optimization procedure every ti oc quantum transitions. In addition, 
if the search looks to be stuck for in a local minimum, we can force a jump 
to another local minimum. 
In Fig. [3]we show the results produced by QSO applied to random instances 

10 



of the Graph Partitioning problem: given a graph G = (V,E), where V de- 
notes the set of vertices and E the set of edges, partition V into two subsets 
such that the subsets have equal size and the number of edges with end- 
points in different subsets is minimized. We extracted our instances from 
the family ^500,0.01 (i- e - random graph with 500 nodes and an edge between 
any two nodes with probability 0.01.) used in Ref.|31| and compare them 
with the results obtained on the same instances by Simulated Annealing. 
The comparison has been made by letting the programs implementing SA 
and QA run for essentially the same (machine) time. Simulated Annealing 
performs better than Quantum Stochastic Optimization: it finds, on the av- 
erage, better (smaller mean value and variance) approximations of the best 
partition. What is interesting is that QSO goes down very fast toward a 
local minimum and then relies on quantum transitions to escape from it. In 
SA, on the contrary, a steep descend toward a local minimum could result 
in a early freezing of the search. We tested QSO also random satisfiable 
instances of 3-SAT with qualitatively similar results. 

Recent variants of QSO, as Imaginary Time Quantum Monte Carlo (ITQMC) 
p6l l45| l46l [38] . introduce a proper annealing schedule in the algorithm: 
the control parameter u(t) is suitably reduced during the algorithm, as 
the temperature is in reduced in SA. Between each "annealing" step, the 
ground state ip u is estimated by the same estimation procedure used in 
QSO. For some optimization problems, such as the Graph partitioning de- 
scribed above, ITQMC performs better than SA, for others, e.g. 3-SAT, 
it does not. There is therefore no a-priori guarantee that QA will produce 
better heuristics than SA [TO] . 

Interestingly enough, it has been shown that ITQMC performs at least as 
well as the real time adiabatic approximation [42j described in the previous 
section. The intuitive reason behind this result has been clearly stated in 
[HUH]: the projection mechanism proper of imaginary time makes the evo- 
lution of the initial condition much more stable than the one determined 
by the Schrodinger equation. The adiabatic theorem, in fact, relates the 
energy gap between the ground state and the first excited state to the an- 
nealing time: the smaller the energy gap, the slower must be the change 
of the Hamiltonian Hj^(t) in order to avoid Landau-Zener transitions of the 
system from the ground to excited states (which do not correspond to the 
solution we are looking for). Since the projection mechanism suppresses 
the amplitude associated to excited states by accidental Landau-Zener tran- 
sitions, the interpolation between the initial Hamiltonian Ha(0) and the 
target Hamiltonian Ht (the annealing) can be carried out more rapidly. 
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Procedure 1 Quantum annealing 



Input: initial condition init; control parameter v\ duration t max ; tunnel 
time t drill] local opt. time ti oc . 

t<-0; 

e <— init] 
Vmin = cost(e); 



repeat 

repeat 

e -(— Quantum Transition(e, v, t max ); 
if cost(e) < v m i n then 

Vmin <- cost(e); 

i,3 <- 0; 
else 

end if 
until i > ti oc 

epsilon -(— Local Optimization(e). 
if cost(e) < v m i n then 
Vmin <- cost(e); 

end if 
until j < tdriii 

draw a trajectory of length vt m ax and jump there. 
Local Optimization(e) 
end while 



Procedure 2 Quantum Transitions 

Input: initial condition e; chain length vt\ set of neighbours to estimate 
Neigh 

for all neighbour k G Neigh do 

estimate the wave function ip u (k); 
end for 

best <— select a neighbour in Neigh with probability proportional to ip v 
return best 
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QSO. Mean = 241.61, StDev= 13.9949 

2d 
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SA. Mean = 226.416, StDev= 1 1.9271 
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Figure 3: Quantum Stochastic Optimization (QSO) vs. Simulated Anneal- 
ing (SA). Benchmark: 100 random instances of ^500,0.01- Allocated CPU- 
time: 600 seconds. Parametrization: QSO: vt = 20; n = 4; \Neigh\ = 5; SA: 
geometric scheduling: temperature at k-th annealing step: t\. = io(0.99) fe ; 
to = 3.0; number of proposed moves per annealing step = 16 • 500. (a) and 
(b) distribution of V m % n for QSO and SA respectively, (c) and (d) Values 
of Vmin reached respectively by QSO and SA on each instance of Graph 
Partitioning after every 50 (machine-time) seconds. In the inset of (c) and 
(d): values of v m i n along a single search path. 

4 Dissipative dynamics 

It becomes quite natural to ask whether a "projection" mechanism similar 
to the one operated by the unphysical imaginary time evolution is available 
in some real time quantum dynamics. If it were possible, we could hope 
to exploit it to speed up quantum adiabatic computation without worrying 
too much about Landau-Zener transitions. Dissipative quantum annealing, 
a novel model proposed in |18|. 116] . represents one first step in answering 
this question. 
Given the usual Hamiltonian H v , we add a non-linear term, the Kostin 
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friction |36|, [37] , modeling the effective interaction of the quantum system 
with an environment which absorbs energy. In the continuous case the 
Schrodinger-Kostin equation reads: 

ij^ = H^ + pK^), (8) 

where H v is the Hamiltonian of ([5]) and 

k<*) = ***{$)■ (9) 

By rewriting the state ip(x) a la de Broglie, 

i;(x) = yW) e lS{x \ (10) 



the nonlinear part K of the Hamiltonian ([Sj) assumes the form K(y/p(x) e w) 
S(x), whose gradient corresponds to the current velocity [H] . This justifies 
the names friction for the term K and friction constant for /3 > 0. 
Given a solution ij){t) of the Schrodinger equation i-iip = (H v + K)ip the 
following inequality holds: 

j t (m\H»\m)<Q- 

It means that the energy of the system is a monotone non increasing func- 
tion of time: the system dissipates energyjj 

In Ref.|18j we showed that friction can play a useful role in suppressing two 
genuinely quantum effects which can affect the exploration of the solution 
space: Bloch oscillations and Anderson localization. Bloch oscillations [12] 
are due to the relation v = sinp between momentum p and velocity v of 
a particle moving on a regular lattice: as the momentum of the particle 
increases monotonically, the velocity can change sing and the particle starts 



4 We point out that the full interaction with the environment would include a ran- 
dom force describing the back action of the environment onto the system. We are then 
well aware of the fact that we are considering just the dissipative part of a fluctuating- 
dissipating system. Indeed, the complete Kostin equation was the correspondent in the 
Schrodinger representation of the Quantum Langevin Equation, that is the equation for 
the observables of an harmonic oscillator coupled to a bath of oscillators [25| in the ther- 
modynamic limit. The proper language to describe such an interaction would be the 
theory of open quantum systems [15] . We decided, however, to adopt a phenomenological 
approach such as the one used in [27] and more recently, in the context of time dependent 
density functional theory, in |52] , A proper analysis of the complete system-environment 
interaction is currently under study. 
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moving backward along the lattice. 

Anderson localization |7J, instead, appears when the lattice presents irreg- 
ularities, manifesting themselves as a random potential on the sites of the 
lattice. The state of a quantum particle moving in a highly irregular poten- 
tial landscape is spatially localized and the probability of tunneling through 
large regions is exponentially suppressed. 

In |18] we considered the simple case of a spin chain of finite size s governed 
by an XY Hamiltonian, that is: 



Hxy = Y, X (^0>i(J + 1) + o-2(j>a(j + !)) 



(11) 



i=i 



We took an initial configuration having a single spin up and all the others 
down. The spin up plays the role of a quantum walker and, by a suitable 
choice of the initial conditions [17], it behaves as an excitation moving ballis- 
tically along the chain as in Figure HI We introduce a simple potential/cost 
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Figure 4: s = 100, A = 1. The excitation, initially localized at the beginning 
of the chain, moves ballistically forth and back. 



function V[x) = —gx. The motion of the walker is affected as shown in Fig- 
Bloch oscillations appear hindering an exhaustive exploration of 
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the solution space. As Figure 5(b) shows, the oscillations can be suppressed 
by adding the frictional term, in the discretized form: 



(K^)(x) = ^Tsm(S(y)-S(y-l)), 



(12) 
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to the system Hamiltonian. Here the phase function S is defined as in 
Eq.dlOD. 

Moreover, the convergence toward the ground state is witnessed by a pro- 
gressive concentration of the probability mass at the end of the chain were 
the ground state is localized. 

Potential profiles determined by optimization problems are usually quite ir- 



regular. Figures 6(a) shows the consequence of the addition of a random 

(ED: 



potential to the free Hamiltonian (|11|) : the probability amplitude remains 
trapped in the first half of the chain, because of Anderson localization. 
Figure |6(b)| shows how dissipation can contrast Anderson localization: the 
probability mass percolates through irregularities and concentrates in the 
region where we chose to localize the ground state by adding a linear poten- 
tial V(x) = —gx as in the previous example. 




4 It' 
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(a) 



(1)) 



Figure 5: s = 100, < t < 20s. Frame (a): g = 6/s, (3 = 0: the wave 
packet gets confined due to Bloch reflection, (b): for g = 6/s, j3 = 8/s: 
Bloch reflections are suppressed by friction. 



5 Conclusion and outlook 

Quantum Annealing is an example of cross contamination between two dif- 
ferent research areas: computer science and physics. The "constructive 
interference" goes in two ways: on one side, a genuinely quantum effect, 
tunneling, suggested a completely classical approximation algorithm which 
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Procedure 3 Local optimization 



Input: initial condition e. 

return the best solution found by any steepest descent strategy. 



40; 
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Figure 6: s = 100, < t < 20s. Frame (a): g = 0, /3 = 0: The wavepacket 
is localized due to Anderson localization (b): for g = 6/s, /3 = 8/s: friction 
suppresses Anderson localization: the probability mass percolates through 
the imperfections (additional random potential extracted from a normal 
population of mean and standard deviation a = 0.06) of the spin chain. 
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provides good heuristics to hard combinatorial problems. In this paper we 
presented in details the earliest version of Quantum Annealing, Quantum 
Stochastic Optimization. 

On the other side, a trick used to estimate the state of minimal energy of 
a quantum system by classical means, Imaginary Time evolution §6§, sug- 
gested the idea of exploiting dissipation to stabilize the evolution of quantum 
systems. 

A mechanism able to suppress Anderson localization would be most welcome 
in Quantum Adiabatic Computation. In fact, Altshuler et al. in Ref.|3j 
advanced the conjecture that "Anderson localization casts clouds over adia- 
batic quantum optimization". The preliminary results presented here on the 
effects Kostin friction on simple toy models are promising: Anderson local- 
ization is suppressed. Still, whether an analogous mechanism can be of some 
use when solving hard optimization problems, together with a quantitative 
assessment of the results presented in [16], is an open research problem. 
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